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ABSTRACT 

We present our cosmic shear analysis of GEMS, one of the largest wide-field surveys ever 
O I undertaken by the Hubble Space Telescope. Imaged with the Advanced Camera for Surveys 

(ACS), GEMS spans 795 square arcmin in the Chandra Deep Field South. We detect weak 
^ ' lensing by large-scale structure in high resolution F606W GEMS data from ~ 60 resolved 

galaxies per square arcminute. We measure the two-point shear correlation function, the top- 
hat shear variance and the shear power spectrum, performing an E/B mode decomposition for 
k> I each statistic. We show that we are not limited by systematic errors and use our results to 

. place joint constraints on the matter density parameter Q„i and the amplitude of the matter 

c3 ■ power spectrum ag. We find = 0.68 ± 0.13 where the la error includes both 

our uncertainty on the median redshift of the survey and sampling variance. 

Removing image and point spread function (PSF) distortions are crucial to all weak lens- 
ing analyses. We therefore include a thorough discussion on the degree of ACS PSF distortion 
and anisotropy which we characterise directly from GEMS data. Consecutively imaged over 
20 days, GEMS data also allows us to investigate PSF instability over time. We find that, even 
in the relatively short GEMS observing period, the ACS PSF ellipticity varies at the level of 
a few percent which we account for with a semi-time dependent PSF model. Our correction 
for the temporal and spatial variability of the PSF is shown to be successful through a series 
of diagnostic tests. 

Key words: cosmology: observations - gravitational lensing - large-scale structure. 



1 INTRODUCTION 

Weak gravitational lensing is a unique probe of the dark mat- 
ter distribution at redshifts z < 1 where, within the cur- 
rently favoured ACDM cosmological paradigm, dark energy be- 
gins to play an important role in the evolution and growth of 
the power spectrum of matter fluctuations. It therefore not only 
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has the power to constrain fundamental cosmological parameters 
such as the matter density para meter ilrn and the amplitude of 
the matter power spectrum ag ("Maoli et al.''200l'; 'Rhode^^lJ 
2001; Van Waerbeke et al. 2001; Hoekstra et al. 2002; Baco n et alJ 
2003; larvisetal. 2003; Brown et al. 2003; Hamana et al. 2003^ 
.Massev et aL2004; Rhodes et aL 2004,; , Van Waerbeke et al,.20oj) . 
but also has the potential to test and constrain quintessence mod- 
els parameterised by the equation of state of the dark energy w{z) 
(Refregier et al. 2004: Jarvis et al]l2005h . Cosmological parame- 
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ter constraints from weak lensing analysis are fully complemen- 
tary to those from cosmic microwave background (CMB) exper- 
iments as parameter degeneracies are almost orthogonal in many 
cases (Brown et al. 2003; Contaldi et al. 2003; Tereno et al. 2005). 
Alone, measurements of the CMB anisotropy at z ~ 1000 are un- 
able to constrain w(z) but combined with future wide-field weak 
lensing surveys, pot entially focused on selected galaxy clusters 
IJain & Tavlodl2003^ ■ the goal of determining w(2:) will certainly 
become attainable. 

Lensing by large-scale structure distorts images of back- 
ground galaxies, inducing weak correlations in the observed ellip- 
ticities of galaxies, termed 'cosmic shear'. The amplitude and an- 
gular dependence of these correlations are related to the non-linear 
matter power spectrum Ps(i) and the geometry of the Universe. 
If we first assume that the minute weak lensing shear distortions 
can be measured in the absence of any systematic errors, arising 
for example from telescope and detector based distortions, the er- 
ror on any weak lensing analysis then has four main sources: shot 
noise from the intrinsic ellipticity distribution of galaxies, sam- 
pling (cosmic) variance, uncertainty in the galaxy redshift distri- 
bution and uncertainty in intrinsic correlations that may exist be- 
tween nearby galaxy pairs. Shot noise can be minimised by survey- 
ing large areas of sky and/or by increasing the number density of 
resolved galaxies in the weak lensing analysis. Sampling variance 
can be minimised by imaging many different lines of sight. With 
spectroscopic or photometric redshift information the galaxy red- 
shift distributio n can be accurately estimated (Brown et al. 2X!Q3.. ; 
IVan Waerbeke e t al. 2005). Redshift in formation also allows for 
the removal of intrinsic correlations teing & Schneideil |2003; 
iHevm ans & Heavens 2003; Hevmans et al. 2004). This effect has 
however been shown to be a small contaminant to the weak lensing 
signal measured from deep surveys (median redshift Zm ~ 1.0), at 
the level of less than a few percent of the shear correlation signal 
teevmans et'd]l2004) . The next generation of weak lensing sur- 
veys aim to obtain precise results, minimising these sources of er- 
ror by surveying of the order of a hundred square degrees with ac- 
companying photometric redshift information. These ground-based 
surveys are however subject to atmospheric seeing which erases the 
weak lensing shear information from all galaxies smaller than the 
size of the seeing disk. This, in effect, limits the maximum depth of 
ground-based weak lensing surveys and hence the sensitivity, lead- 
ing to proposals for future deep wide-field space-based observa- 
tions (Rhodes et al. 2004). With space-based data the number den- 
sity of resolved galaxies, in comparison to ground-based studies, 
is multiplied two-fold and more. Resolving images of more dis- 
tant galaxie s will permit high r esolution maps of the dark matter 
distribution ^ Massev et all2004l ) and will yield significantly higher 
signal-to-noise constraints on cosmological parameters in compar- 
ison to constraints obtained from the same survey area of a seeing 
limited ground-based survey (Brown et al. in prep). 

With the installation of the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope (HST), relatively wide-field 
space-based weak lensing studies are feasible and in this paper 
we present the detection of weak gravitational lensing by large- 
scale structure in GEMS; the Galaxy Evolut ion from Morph ology 
and Spectral energy distributions survey (see lRix et all2004l for an 
overview). Spanning 795 square arcmin, GEMS is currently the 
largest space-based colour mosaic. With a high number density of 
resolved sources and spectroscopic and/or photometric r edshifts for 
^ 80 00 of our sourc es from the COMBO - 17 survey jWolfetalJ 
|2004 and the VVDS <Le Fevre et alj2004 . we can beat down shot 
noise and obtain a good estimate of the galaxy redshift distribution. 



The error sources discussed above have not included system- 
atic errors, but in reality for all weak lensing surveys the accuracy 
of any analysis depends critically on the correction for instrumen- 
tal distortions, which are orders of magnitude larger than the un- 
derlying gravitational shear distortions that we wish to detect. The 
strongest distortion results from a convolution of the image with 
the point spread function (PSF) of the telescope and camera. For 
the detection of cosmic shear at the level of only a few percent, the 
PSF distortion needs to be understood and controlled to an accu- 
rac y of better than one percent. For space-based weak lensing stud- 
ies ('Hoekstra et al.'l998^.'Rhodes et al.'200C|.|200l|;|Refregier et alJ 
^002; Hammerle et al. 2002; Casertano et alj |2003t iRhodes et 4] 
2004; Miralles et al. 2005), the limited number of stars in each 
camera field of view provides insufficient coverage to map the PSF 
anisotropy. This has made it necessary to assume long-term sta- 
bility in the PSF pattern, often using PSF models derived from 
archived images of globular clusters, with low order focus correc- 
tions, or from Tiny Tim (Rrist 2000) thereby critically limiting the 
accuracy of the correction for PSF distortions. In contrast to previ- 
ous space-based lensing studies we are not forced to assume long- 
term PSF stability as the GEMS PSF can be characterised directly 
from the wide-area data where all but three out of the sixty-three 
ACS images were observed in the space of 20 days. Hence we need 
only assume short term PSF stability, which we test in Section l4m 
A detailed investigation into the GEMS ACS PSF will be presented 
in Jahnke et al. (in prep) with a view to many different astronomical 
applications. 

This paper is organised as follows. In Section|2|we summarise 
the basic theory that underpins cosmic shear studies and review 
the particulars of our analysis. We discuss, in Section |3| the ob- 
servations and data reduction, paying specific attention to potential 
sources of image distortions. We focus on characterising and cor- 
recting for the temporally and spatially varying dominant PSF dis- 
tortion in Section|4]and compare, in Section|5| galaxy shear mea- 
sured from the two GEMS passbands; F606W and F850LP. In Sec- 
tion |6| we estimate the redshift distribution of the GEMS survey 
from COMBO- 17 and VVDS redshift catalogues. Our statistical 
analysis is detailed in Section |7| where we measure the shear cor- 
relation function, the top-hat shear variance statistic and the shear 
power spectrum performing many diagnostic tests for systematic 
errors. These results are used to place joint constraints on Q.m and 
(Tg in Section|8| We conclude in Section|9|with a comparison to re- 
sults from other space-based and ground-based cosmic shear anal- 
yses. 



2 METHOD 

The theory of weak gr avitational le nsing by large structure, detailed 
in Bartelmann & Sch neideJ i200lh . directly relates the non-linear 
matter power spectrum Ps to the observable weak lensing (com- 
plex) shear field 7 = 71 + 272 characterised through the shear 
(convergence) power spectrum Pk, 

P.„.5^f2kr.„s!Mp, (''„■). ,„ 

4c4 a^{w) \fK(w) J 

where w is the comoving radial distance, a{'w) is the dimension- 
less scale factor. Ho is the Hubble parameter, Qrn the matter den- 
sity parameter and g{w) is a weighting function locating the lensed 
sources, 

g(w)= / dw (l>[w) — ; S (2) 
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where (f)[w{z)] dw is the o bserved number of ga laxies in dw and 
wh is the horizon distance iSchneider et alll998h . Note in this pa- 
per we will assume a flat Universe with zero curvature where the 
comoving angular diameter distance = w. 

In this paper we measure directly the shear power spectrum 
Pk, the shear correlation function {"f {0)^(0 + AO)) which we split 
into a tangential component 



<7t7t)e = 4^ y deiPM) [Jo{ee) + MiO)], 
and a rotated component 
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and the top-hat shear variance {\j\'^)g measured in a circle of an- 
gular radius 



{7^)9 = 
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(5) 



Note Ji is the i* order Bessel function of the first kind. 

In order to exploit the straightforward physics of weak lens- 
ing, one requires an estimate of the gravitational shear experienced 
by each galaxy . Kaiser et al. ( 1995), LuDoino & Kaiser ( 1997) and 
iHoekstra e t al. ( 1998) (KSB-I-) prescribe a method to invert the ef- 
fects of the PSF smearing and shearing, recovering an unbiased 
shear estimator uncontaminated by the systematic distortion of 
the PSF. Objects are parameterised according to their weighted 
quadrupole moments 



J d^6W{e)i(p)eiej 
J d'^ew{e)i{e) 



(6) 



where / is the surface brightness of the object, 9 is the angular dis- 
tance from the object centre and is a Gaussian weight function 
of scale length rg, where Vg is some measurement of galaxy size, 
for example the half light radius. For a perfect ellipse, the weighted 
quadrupole moments are related to the weighted ellipticity param- 
eters by 
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lKaiseretalJ<1995h show that if the PSF distortion can be described 
as a small but highly anisotropic distortion convolved with a large 
circularly symmetric seeing disk, then the ellipticity of a PSF cor- 
rected galaxy is given by 



cor obs 



(8) 



where p is a vector that measures the PS F anisotropy, and P"™ 
is the smear polarisability tensor given in 'Hoekstra et alj jl998l) . 
p{9) can be estimated from images of stellar objects at position 
by noting that a star, denoted throughout this paper with *, imaged 
in the absence of PSF distortions has zero ellipticity: 6'^'^°'^ — 0. 
Hence, 



/ 7->sm* \ — 1 * 



(9) 



For space-based imaging, where PSFs deviate strongly from a 
Gaussian, this PSF correction is mathematically poorly defined 
^Kaiser 2000) such that it is important to calculate the PSF correc- 
tion vector p not only as a func tion of galaxy p osition 0, but also 
as a function of galaxy size rg ^Hoekstra et al. 1998;). This rather 
unsatisfactory situation has prom pted the develop ment of altema- 
tive methods jRhodes et all200(]tlKaiseil200a:lBernstein & JarvisI 



l2002tlRefregier & Bacoi]l2003l : lMassev & Refregiej2004 but for 
the purpose of this paper we will focus on the most commonly used 
KSB+ technique, deferring our analysis with different techniques to 
a future paper. 

The isotropic effect of the PSF is to convolve galaxy images 
with a circular kernel. This makes objects appear rounder, erasing 
shear information for galaxies smaller than the kernel size, which 
have to be removed from the galaxy shear sample. For the larger 
galaxies, this resolution effect can be accounted for by applying the 
pre-seeing shear polarisab ility tensor correction P'' , as proposed by 
iLuppino & Kaised such that 



cor s 1 T-j'y 



(10) 



where e" is th e true source ellipticity an d 7 is the pre-seeing gravi- 
tational shear. iLuppino & Kaised iT997ft show that 
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where P'^^ is the shear polarisability tensor given in lHoekstra et alJ 
1*1998^ and P""* and P'''" are the stellar smear and shear polaris- 
ability tensors respectively. This relation is only strictly true when 
all values are measured from the PSF de-convolved image which 
is difficult to create in practice. P^ is therefore calculated from the 
PSF distorted images which produces very noisy measurements. 

Combining the PSF correction, equation Js}, and the P^ see- 
ing correction, the final KSB-F shear estimator 7 is given by 



(12) 



When averaging over many galaxies, assuming a random distribu- 
tion of intrinsic galaxy ellipticities, (e^) = 0, and hence (7) = 7, 
providing a good estimate for the gravitational shear. 



3 THE GEMS DATA 

The GEMS survey jRix et alj20'04l) spans an area of ~ 28' x 28' 
centred on the Chandra Deep Field South (CDFS), combining 125 
orbits of ACS/H ST time with su pplementary data from the GOODS 
project ( Giavali sco et all2004b . 78 ACS tiles have been imaged in 
F606W and 77 ACS tiles in F850LP, where the point source 5a 
detection limits reach meoe = 28.3 and mgso = 27.1. In this sec- 
tion we review the data set, discuss the potential for bias within the 
source catalogues and highlight possible sources of image distor- 
tion, in addition to the strong anisotropic PSF distortion which is 
characterised and corrected for in Section|4] A detailed account of 
the full GEMS data reduction method will be presented by Cald- 
well et al. (in prep). 

The ACS wide-field camera has a field-of-view ~ 3.4 x 3.4 
arcmin com prising two 4 096 x 2048 CCD chips of pixel scale 
0.05 arcsec iFord et aL 2003). GEMS observes, in sequence, three 
separate exposures per ACS tile dithered by ~ 3 arcsec, where 
the observation strategy has been designed in such a way so as to 
bridge the inter-chip gap and provide sub-pixel oversampling of 
pixel scale 0.03 arcsec in the final co-added image (~ 7000 x 7000 
pixels). GOODS have employed a different observing strategy, us- 
ing only two separate dithered exposures. In ord er to optimise the 
survey for Supernova searches iRiess et all2004> the GOODS area 
is re-imaged in 5 different epochs, but to obtain similar depths 
to the GEMS data and minimise the effects of PSF time varia- 
tion, we co-add only the exposures from the first epoch of ob- 
servations. This however leaves us with a slightly shallower cen- 
tral region in the F606W GEMS mosaic which can be seen from 
the median magnitude of each data set; m606(GEMS) — 25.6, 
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meoe (GOODS) = 25.1. We note that the 2 exposure GOODS 
dithering pattern will result in a poorer cosmic ray rejection in the 
GOODS area of the GEMS mosaic, and will impact somewhat on 
the PSF. 

Images from the ACS suffer from strong geometric distortions 
as a result of the off-axis location of the camera within HST. In ad- 
dition, the HST optical assemb ly and the ACS mirrors also induce 
distortions. Meurer et al. ( 2003) accurately calibrate and model this 
distortion from dithered images of star clusters. With this model all 
tiles are drizzled onto a celestial pixel grid us ing a version of the 
multidrizzle software iKoekemoer et alJl2003h . where the astrom- 
etry of each GEMS tile is tied to the ove rall catalogue from the 
ground-based COMBO-17 R band image fWolf et alj |2'004h . Av- 
eraging our final PSF corrected shear catalogues as a function of 
ACS chip position, we find no evidence for any residual geometric 
distortions remaining in our final multi-drizzled images. 

A second order geometric distortion arises from the effect of 
velocity aberration l lCox & Gi lliland 2002). The HST guides on 
nearby stars peripheral to the ACS field of view, making small cor- 
rections to keep the primary target on a fixed pixel. It cannot how- 
ever correct for the isotropic plate-scale breathing which is of the 
order ~ 0.001% over an orbit. These small changes in pixel scale 
are allowed for and corrected by our version of multidrizzle, pro- 
vided the individual exposures are on the same scale. Variation of 
scale during each observation would result in a slight blurring of 
the co-added images that increases radially from the centre. GEMS 
ACS tiles, observed over 1/7 of the orbit, suffer from pixel scale 
variation that is, at maximum, a difference of 0.004 arcsec corner 
to comer. Our results later show no significant variation in aver- 
age galaxy shear as a function of chip position and we hence con- 
clude that this effect is not significant within our measurement ac- 
curacy. In future data reductions of GEMS we will include a veloc- 
ity aberration correction to the pixel scale using an updated version 
of multi-drizzle. 

We use the SExtractor software iSertin & Arnoutslll996h to 
detect sources on both the F606W and F850LP imaging data, with 
the two-step approach described in Rix et al.. (.2004,) . In short, this 
method of combining two SExtractor catalogues, one measured 
with high signal-to-noise detection thresholding, and one measured 
with low signal-to-noise detection thresholding, allows us to find 
the best compromise between detecting faint galaxies without de- 
blending bright nearby galaxies into many different components. 
We also use SExtractor to determine the weak variation of the sky 
background across the tile, which is then subtracted for the follow- 
ing KSB-l- analysis. We define galaxy size Vg as the half light radius 
measured by SExtractor (f lux.radius)^ and calculate weighted 
ellipticity parameters Si and the shear and smear polarisability ten- 
sors; P^^ and P"™ for each object in the SExtractor catalogue. 

The accuracy of the centroid determined by SExtractor 



^ SExtractor analysis of image simulations have shown that the 
flux.radius does not accurately measure the true half-light radius 
(HauBler in prep.). A good relationship between f lux.radius and the 
input half light radius is seen but this relationship differs between disk and 
bulge dominated galaxies. The value that we choose for Tg does not change 
the results as the tinal P"' correction, equation 1111 . includes a correction 
for the Tg Gaussian weighting. The choice of rg does however affect the 
noise on each measurement and so we choose Vg = /x f lux.radius 
where / is chosen to minimise the noise in the measured galaxy ellipticities. 
We find / = 1. Note that alternative ways of defining rg using SExtractor 
ellipticity measures (see for example Rhodes et al. 2000) have been found 
to yield slightly noisier results with the GEMS data. 



is directly linked to the accuracy of each KSB-l- galaxy shear 
measurement. In the presence of non-isotropic c entroiding er- 
rors there is also the potential for centroid bias (Kaiser 2000; 
Bernstein & Jarvis 2002) which can arise if, for example, the 
errors in the x direction exceed those in the y direction producing 
a tendency to bias towards an average galaxy shear in the x 
direction. The GEMS galaxies have been modelled usi ng the two- 
dimensional galaxy profile model fitting code GALFIT ('Pe ng et all 
.200Z) which finds the best fit PSF convolved Sersic profiles to 
eac h galaxy, allowing th e centroid to be a free parameter in the fit 
("see lBarden et alj|2004l and HauBler et al. in prep for details). We 
can therefore test if we are subject to centroid bias by comparing 
GALFIT and SExtractor centroids. The average pixel offset is 
consistent with zero in the x direction, Aa; = —0.001 ± 0.02 
and very close to zero in the y direction Ay = 0.03 ± 0.02. 
Calculating the ellipticity of a mock circular Gaussian galaxy 
A'^ — 10* times, assuming Gaussian distributed centroid errors 
with mean and width as estimated by GALFIT, we find that for our 
smallest galaxies, centroid errors induce a systematic centroid bias 
[ei = (-2.754 ±0.001) x 10"*, 62 = (-7.14 ±0.02) x lO"*^], 
which is negligible compared to our current measurement accuracy. 
For larger galaxies centroid bias decreases. Note that the GEMS 
GALFIT galaxy profile parameters cannot currently be used for 
weak lensing studies as the PSF has been derived from co-added 
stellar images and therefore does not allow for the spatial variation 
of the distortion. With a spatially varying PSF model GALFIT 
could be used for measuring galaxy shear, although this would be 
a time consuming process. 

When compiling source catalogues one should consider selec- 
tion bias where any prefer ence to se lect galaxies oriented in the 
same direction as the PSF (Kaiser 2000) and galaxies that are anti- 
correlate d with the gravitationa l shear (and as a result appear more 
circular) iHirata & Seliakl2003h . would bias the mean ellipticity of 
the population. Through simulations of artificial disc galaxy light 
profiles, convolved with the ACS PSF (see HauBler et al. in prep 
for details), we see no significant selection bias when we intro- 
duce a SNR > 15 selection criteria (defining SNR = flux / 
f lux.error), i.e within the noise of the sample, there is neither 
a preference for selecting faint galaxies oriented with the PSF, nor 
a preference for selecting more circular faint objects. 

To remove erroneous detections along the chip boundaries, 
diffraction spikes from stars, satellite trails and reflection ghosts, 
each imag e catalogue is masked b y hand using the method de- 
scribed in MacDonald et alJ J2004h . Using weight maps to define 
the best regions in each tile, we combine the masked catalogues 
from each ACS tile ensuring that in the overlapping regions of 
neighbouring tiles, only the data from the better tile is included. 
Note that the objects in overlapping regions are used for consis- 
tency checks to test the accuracy of the galaxy shear measurement. 

The charge transfer efficiency (CTE) of space-based instru- 
ments provides another source of image distortion. Objects with 
low signal-to-noise in low sky background images tend to bleed 
in the read-out direction of the CCD camera, causing an elonga- 
tion of the objects that is correlated with the read-out direction and 
the distance from the read-out amplifier. Over time the CTE de- 
grades, increasing t he magnitude of this effect ( Mirall es et alj20o3: 
iRhodes et all2004^ . The read-out-amplifiers for the ACS lie at each 
comer of the camera, with the read-out direction along the y axis 
of the CCD and we find no correlation between the average PSF 
corrected galaxy shear along the y axis, and the galaxy distance 
from the read-out-amplifiers. Hence, even though the CTE of the 
ACS wide-field camera has been shown to degrade with time iRiessI 
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|2002» . we find no signature for CTE in our data, potentially a result 
of the data being observed relatively soon after the ACS installa- 
tion. 



4 CHARACTERISATION OF, AND CORRECTION FOR, 
THE ANISOTROPIC ACS PSF 

In this section we give a thorough account of the techniques which 
we have used to characterise and correct for the anisotropic ACS 
PSF. To date there has not been such a large set of HST data im- 
aged in a short time frame which can allow for a rigorous semi-time 
dependent PSF analysis. The PSF is characterised through images 
of non-saturated stellar objects that are selected through their lo- 
cus in the size-magnitude plane. Stars can be easily identified, ex- 
tending from bright magnitudes to faint magnitudes into the main 
distribution of galaxies, remaining at one characteristic size. We 
use both the full width half maximum (FWHM) and the half light 
radius (f lux_radius) measured by SExtractor as a definition 
of size to select the stellar objects, where stellar candidates must 
lie along the stellar locus in both the FWHM-magnitude plane and 
the f lux.radius-magnitude plane. This method selects ~ 1300 
stellar objects in the F850LP images and ~ 900 stellar objects 
in the F606W images. Note that the loss of stellar objects in the 
F606W arises from the increased number of saturated stars and the 
use of more conservative cuts at fainter magnitudes to avoid confu- 
sion with small faint galaxies. 

For each star we measure the weighted stellar ellipticity pa- 
rameters Ea and the stellar smear polarisability tensor p=™* using 
Gaussian weights W{rg) with different smoothing scales; rg rang- 
ing from 1.9 pixels, the minimum stellar size measured by SEx- 
tractor, to 10 pixels. We limit the maximum smoothing scale to 
avoid excessive inclusion of light from neighbouring objects which 
quickly introduces noise into the stellar shape measurement. Fig- 
ure Q shows the variation in the stellar ellipticity parameters 
across the ACS field of view, measured using all the GEMS F606W 
data, (upper panels) and all the GEMS F850LP data (lower pan- 
els). The horizontal spacing at y ^ 3200 pixels results from the 
chip boundaries where shape estimates become unreliable. For Fig- 
ure Q £* has been measured using two differently scaled weight 
functions W{rg), the first looking at the core PSF distortion with 
Vg = 2.5 pixels (left panels), and the second looking at the PSF dis- 
tortion averaged over the main extent of the star, with rg = 7 pixels 
(right panels). This figure shows that the PSF distortion is clearly 
anisotropic and varies with scale size, and with filter. The F850LP 
PSF has a strong horizontal diffraction spike, which dominates the 
average PSF distortion on large scales resulting in a strong positive 
£i component across the full field of v iew. This d i ffracti on spike 
may therefore account for the claim bv lPark et alj i2004h that the 
ACS PSF is fairly isotropic. 

In order to accurately characterise the anisotropy of the PSF 
across the field of view of each ACS tile in the survey we wish 
to maximise the surface density of stellar objects as a function of 
{x, y) position. This however necessitates some assumptions about 
the PSF stability over time as the selected stellar number density 
corresponds to only ~ 16 stars per ACS tile in F850LP and ~ 11 
stars per ACS tile in F606W. Figure Q shows smooth variation in 
the PSF as a function of chip position indicating that any varia- 
tion of the PSF in the 20 day duration of the GEMS observations 
is small. We split our stellar sample into stars imaged by GEMS 
and stars imaged by GOODS, as each data set derives from co- 
added exposures with different dithering patterns which impact on 
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Figure 1. The anisotropic ACS PSF measured from stellar sources in 
F606W images (upper panels) and F850LP images (lower panels). The stel- 
lar ellipticity e* plotted has been measured using two differently scaled 
weight functions Wivg); left panels Vg = 2.5 pixels (PSF core distortion), 
right panels Tg = 1 pixels (PSF wing distortion). The 5% bar in the upper 
left corner of each panel shows the scale, which is the same for each panel. 



the PSF. The first epoch of GOODS observations spanned 5 days, 
and all but 3/63 GEMS tiles were observed in the space of 20 days. 
We reject from our analysis the 3 GEMS tiles which were taken 
out of sequence and split our GEMS sample into 2 data sets as- 
suming PSF stability on the scale of 10 days. We will quantify the 
validity of this assumption in Section liHI Our PSF model is thus 
semi-time dependent as we use three time-dependent PSF models 
derived from fits to three different data sets selected by observation 
date. 

To model the anisotropy of the PSF across the field of view, 
we fit a two-dimensional second order polynomial to the PSF cor- 
rection vector p equation J9}, modelling each CCD chip and data 
set separately. Before fitting we remove outliers with a 3cr devia- 
tion from py^, and then iterate twice during the fit, removing out- 
liers with a 3(T deviation in their PSF corrected ellipticity from 
ejj""^*. Figure |2| shows the variation in the measured PSF correc- 
tion vector p across the ACS field of view, measured using all 
the GEMS F606W data (upper left). The p values calculated from 
our semi-time-dependent polynomial models (upper right), the cor- 
rected stellar ellipticities e™"^* (lower left, note that p and e are 
plotted on different scales), and the ellipticity distribution of stars 
before and after the PSF correction (lower right) are also shown. 
For this figure we have used a smoothing scale of Tg = 5.9 pixels, 
which is the median galaxy size in our catalogue. In this case we 
find the mean stellar ellipticity before and after correction to be 

e\ = 0.0475 ± 0.015 el = 0.0157 ± 0.014 (before), 
e\ = 0.0003 ± 0.0007 el = -0.0001 ± 0.0007 (after). (13) 

This demonstrates that the PSF correction significantly reduces the 
mean stellar ellipticity such that it is consistent with zero, and that 
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Figure 2. The upper left panel shows the variation in the measured PSF 
correction vector p across the ACS field of view, measured using all the 
GEMS F606W data. This data is modelled with a semi-time dependent two- 
dimensional second order polynomial shown in the upper right panel. The 
PSF corrected stellar ellipticities e™"^*, lower left panel, display random 
orientation. The 5% bar in the upper left corner shows the scale. Note that 
p and e are not directly equivalent and are thus plotted on different scales. 
The ellipticity distribution of stars before (circles) and after (dots) PSF cor- 
rection are also shown (lower right). 



the dispersion also decreases by a factor ~ 2. Note that the success 
of this correction does lessen somewhat with increasing rg, as the 
noise in the measurement of p grows, but as the number of galax- 
ies to which these high Vg corrections apply decreases in turn, this 
effect is not problematic. 

With the PSF models we correct our galaxy catalogue for PSF 
distortions using equation js}. To test for residual PSF related sys- 
tematic distortions we search for correlations between PSF cor- 
rected galaxy ellipticity e™"^ and stellar ellipticity e* . Our first sim- 
ple test splits the survey into square cells of side 20 arcseconds 
in the ACS {x, y) plane, calculating the cell averaged PSF cor- 
rected galaxy ellipticity {ef"^) and the cell averaged uncorrected 
stellar ellipticity (e*) determined at rg = 5.9, the median galaxy 
size in the survey. Figure|3|(lower panel) shows the resulting mean 
Tjl^ef^) /N cells as a function of cell stellar ellipticity (e*), where 
for comparison we also show the average ellipticity of galaxies, 
which have not been PSF corrected, as a function of stellar ellip- 
ticity (upper panel). The correlation found with the uncorrected 
galaxy catalogue is not seen in the PSF corrected galaxy elliptic- 
ities, indicating the success of the PSF correction. This promising 
result will be tested more rigorously in Section2j2| 

4.1 The temporal stability of the PSF 

PSF time variation in space-based instruments is known to result 
from telescope 'breathing', as the HST goes into and out of sun- 
light in its 90 minute orbit, a nd from a slow chan ge in focus which 
is periodically corrected for iRhodes et al]|200(]l) . Variation in the 
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Figure 3. The upper panel shows the correlation between the mean ob- 
served galaxy ellipticity (e°'"') averaged in square cells of side 20 arcsec- 
onds, and the cell stellar ellipticity (e*) measured with rg = 5.9 pixels. 
CoiTection for the anisotropic PSF removes this congelation which can be 
seen in the lower panel where the mean PSF corrected galaxy ellipticity 
(e^""'} is shown as a function of stellar ellipticity. Note that the upper and 
lower panels are plotted on different scales. 



PSF as measured from reduced images can also be caused by slight 
differences in the data reduction method but the consistent GEMS 
observation and reduction strategy minimises this effect. With our 
large set of HST data we are able to test the stability of the ACS PSF 
by looking for changes in the average stellar ellipticity as a func- 
tion of time and time dependent changes in the anisotropy of the 
PSF distortion. Figure|4]shows the variation in the average F606W 
and F850LP stellar ellipticity parameters as a function of observa- 
tion date^ where for each ACS image the F850LP and F606W data 
were taken in succession. This figure reveals a clear trend in both 
filters with el (circles) increasing and £2 (squares) decreasing by a 
few percent during the observation period, a variation that is of the 
order of the signal we wish to detect. This measurem ent of ACS 
PSF temporal instability is in agreement with ljee et aL (2005) who 
show that their PSF can only be characterised from archived stellar 
cluster images when a small ellipticity adjustment is applied. 

The reason for the temporal variation of the ACS PSF is not 
fully understood. Figure |4] shows that the F606W PSF becomes 
more circular (|e| ^ 0) in time, in contrast to the F850LP PSF 
which becomes more elliptical. This could potentially be explained 
by a slow de-focus if one also considers the poorly understood 
strong horizontal diffraction spike seen in the F850LP PSF. It is 
unlikely that this spike is caused by the ACS optics and it will there- 
fore remain unaffected by any de-focus. A de-focus will circularise 
the images, as seen from the F606W data, and lower the contrast 
between the PSF core and diffraction spike in the F850LP data, 
increasing the F850LP el component. 



^ Our consistent GEMS reduction strategy minimises the chances that the 
PSF variation we see results from data manipulation but the reader should 
note that the observation strategy of GEMS does correlate observation date 
and declination. 
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Figure 4. The stellar ellipticity parameters (circles) and (squares) av- 
eraged across the full ACS field of view as a function of observation date. 
The stellar ellipticity averaged over the whole data set is plotted as an iso- 
lated point on the far left of the plot. The upper panel shows the variation 
measured from F6()6W data and the lower panel shows the variation mea- 
sured from F85()LP showing that both filters follow the same trends in the 
and £2 components. The eiTors plotted are the errors on the mean stel- 
lar ellipticity. For this figure we have used a smoothing scale of Vg = 5.9 
pixels to calculate the stellar ellipticities, but we see PSF variation over the 
full range of smoothing scales that we use in this analysis. 



Figure |2| shows that our semi-time-dependent PSF correction 
reduces the average stellar ellipticity to zero and this generally 
holds when we measure the PSF corrected stellar ellipticity as a 
function of observation date, as shown in Figure|5|- There are how- 
ever notable exceptions, for example ACS images taken on the third 
and sixth observation days, showing that the semi-time-dependent 
modelling is only a partial solution for time varying PSF modelling. 

To investigate the effect on the anisotropy of the PSF Figure|6| 
compares the difference between the stellar ellipticity predicted by 
the PSF models calculated for the two GEMS data sets (days 1-10 
and days 11-20). To convert Ap into a difference in stellar ellip- 
ticity we multiply by the average stellar smear polarisability p=™* 
following equation (S). This comparison reveals variation in the 
anisotropy of the PSF at the level of, at maximum, Ae* ^ 5%. We 
hence conclude that the PSF time variation is not simply a change 
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Figure S. The F606W corrected stellar ellipticity parameters e (circles) 
and e™'* (squares) averaged across the full ACS field of view as a function 
of observation date. The coirected stellar ellipticity averages to zero over 
the whole data set as shown by the isolated points on the far left of the plot. 
In general the PSF correction works well reducing the corrected ellipticity 
as a function of observation date close to zero, and significantly reducing the 
scatter (i.e removing the anisotropy) seen in Figure|3(plotted on a different 
scale). There is however the occasional case where the stellar ellipticity 
is significantly different from zero showing that the semi-time-dependent 
modelling is only a partial solution for time varying PSF modelling. The 
eiTors plotted are the errors on the mean corrected stellar ellipticity. 



in average ellipticity but rather an instability in the PSF anisotropy 
across the field of view. 

Creating PSF models for each ACS tile based on at maximum 
a few tens of stellar objects will yield systematic errors larger than 
the variation in the mean corrected stellar ellipticity that we see 
in Figure |5| Thus, the short-term PSF variation cannot be simply 
modelled with the method that we have used. The time dependent 
variation of the PSF whilst significant is at a low level and we there- 
fore proceed with our semi-time dependent PSF model noting that 
short term variation in the PSF may well contribute to systematic 
errors (see Section l731 . 



4.2 Isotropic correction 

The application of the anisotropic PSF correction to observed 
galaxy ellipticities, through equation Js}, leaves an effectively 
isotropic distortion. This distortion makes objects rounder as a re- 
sult of both the PSF and the Gaussian weight function used to mea- 
sure the galaxy shapes, affecting smaller galaxies more strongly 
than the larger galaxies. To correct for this effect and convert 
weighted galaxy ellipticities e into unbiased shear estimators 7, 
we use the pre-seeing shear polarisability tensor P^, equation 
<1 H . is calculated for each galaxy from the measured galaxy 
smear and shear polarisability tensors, P°™ and P^^, and a term 
that is dependent on stellar smear and shear polarisability tensors; 
^psmt-j-i p|h*^ jjjg p7 correction is isotropic we can calcu- 
late this stellar term purely as a function of smoothing scale rg, 
averaging over all the stellar objects that were used in the previous 
anisotropic PSF modelling. Note that we calculate a different value 
for each data set as the PSF variation we see may well be related 
to camera focus which will effect the P^ correction as well as the 
anisotropic PSF correction. P^ is a very noisy quantity, especially 
for small galaxies, but as we expect there to be no difference in the 
P^ correction for the 71 and 72 shear components, we can reduce 
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Figure 6. The difference between stellar ellipticity calculated from PSF 
models for each of the two GEMS data sets with Vg = 5.9 pixels. This 
shows that PSF variation is not a simple change in average ellipticity but 
also an instability in the PSF anisotropy across the field of view, varying at 
the maximum level of Ae* = 5%. 
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Figure 7. Comparison of galaxy shear 7i measured in both F606W and 
F850LP GEMS imaging. The grey-scale shows the number density of ob- 
jects which cluster at low shear values (black = 210 galaxies, white — > 
galaxies). Over-plotted are the grey-scale contours which follow the 1:1 re- 
lationship that we would wish to see between the two data sets. 



this noise somewhat by treating P'' as a scalar equal to half its trace 
(note that the off diagonal terms of are typically very small). 

In an effort to reduc e the noise on still further, is of- 
ten fit as a function of jHoekstra et alll998tlBacon et alllOol 
iBrown et all2003HMassev et all2004) . With space-based data this 
fitting method produces a bias as P'' is a strong function of galaxy 
ellipticity where the dependence can be demonstrated by consid- 
ering that in the case of a galaxy observed in the absence of PSF 
smearing and shearing, P^ reduces to P"' — 2(1 — e^), where 
e is the unweighted galaxy ellipticity. 7 — e/P^ is very sensi- 
tive to small errors in a functional fit of P~'{rg,e) and we there- 
fore do not use any form of fitting to P'' . Although this decreases 
the signal-to-noise of the shear measurement, it avoids any form of 
shear calibration bias which would not be identified with an E/B 
mode decomposition as discussed in Sections [7.2.2l and l7.4l 



4.3 Catalogue selection criteria 

For our final PSF corrected shear catalogue we select galaxies 
with size rg > 2.4 pixels, galaxy shear I7I < 1, magnitude 
21 < mgoe < 27, and SNR > 15. We remove galaxies from 
the catalogue with neighbouring objects closer than 20 pixels (0.6 
arcsec) to reduce noise in the ellipticity measurement from over- 
lapping isophotes. These selection criteria yield 47373 galaxies in 
the F606W images, a number density of ~ 60 galaxies per square 
arcminute, and 23860 galaxies in the F850LP images and a num- 
ber density of ^ 30 galaxies per square arcminute. Note ^ 15% 
of our F606W sources have photometric redshift estimates from 
COMBO- 17. We find no significant correlations of the galaxy shear 
with chip position, galaxy size, magnitude or SNR. 

To analyse the full GEMS mosaic we rotate the shear mea- 
surements from each ACS tile into a common reference frame by 

( iT' \^( cos 20 sin2<^ \( n4^ 
72°' ) \ -sm2^ cos2,^ ) \ -^2 J' ^ ' 



where (j> is defined to be the angle between the x axis of each ACS 
tile and a line of constant declination. 



5 COMPARISON OF r606W AND F850LP DATA 

As the F850LP data is significantly shallower than the F606W data 
we omit it from our cosmic shear analysis, but it is interesting to use 
the galaxies detected in both F606W and F850LP as a consistency 
check to test if our method is sensitive to the differences in filter 
PSFs seen in FigureQ Even though we expect galaxy morphology 
to appear differently in the F606W and F850LP images, we expect 
our shear estimates to remain consistent. Figure |7| shows the very 
good agreement between galaxy shear measured in the F606W and 
F850LP images where the grey-scale shows the number density of 
objects. This comparison shows that our method of measuring shear 
produces very consistent results for galaxies imaged with different 
PSFs and different noise properties, showing no significant calibra- 
tion biases. In principle one should correct the galaxy ellipticity 
based on galaxy colour because of the chromatic anisotropy of the 
PSF, but this comparison also shows that the colour of the PSF does 
not significantly impact on the shear measurement. 



6 GEMS REDSHIFT DISTRIBUTION 

To interpret a weak lensing signal we need to know the redshift 
distribution of the lensed sources (see equation J2j). The deeper 
the survey is, the stronger the signal we expect to measure. We 
estimate the source redshift distribution of GEMS bas ed on pho- 
tomet ric redshifts from the CDFS COMBO-I7 survey jWolf et alj 
|2004|) and spec troscopic redshifts from the CDFS VVDS survey 
iLe Fevre et ai] Eo04). by assuming that a magnitude dependent 
redshift distribution can be parameterised as 
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Figure 8. The median galaxy redshift as a function of msoe magnitude 
based on photometric redshifts from COMBO-17 (squares) and the HDFN 
(triangles), and spectroscopic redshifts from VVDS (circles). Median red- 
shift errors include uncertainty in redshift completeness which dominates 
the COMBO-17 and VVDS results at faint magnitudes, and statistical un- 
certainty which dominates in the smaller galaxy samples at bright magni- 
tudes. For the photometric redshifts we also include the average redshift er- 
ror, added in quadrature. Magnitude errors show the uncertainty in the mean 
msoe magnitude in each magnitude bin selected in Iab (VVDS), Rvoga 
(COMBO-17) and meoe (HDFN). Over-plotted is the best linear fit to the 
COMBO-17 and VVDS data. The HDFN data is shown at faint magnitudes 
to justify our extrapolation to faint magnitudes of the COMBO-nA'VDS 
linear fit. 



n(z, mag) oc 2 exp 



2:0 (mag) 



(15) 



where zo is calculated from the median redshift Zm with zq = 
2^/1.412 (Bauah & Efstathiou 1994). We calculate 2,„(mag) for 
each survey in magnitude bins of width 0.5 magnitude, out to the 
limiting magnitude of each survey VVDS, Iab < 24; COMBO-17, 
Rvega < 24. These estimates are taken as lower limits for the true 
median redshift, as both surveys suffer from redshift incomplete- 
ness at faint magnitude s. To calculate uppe r limits for the true me- 
dian redshift we follow Isrown et all assuming those galax- 
ies without an assigned redshift are most likely to be at a higher red- 
shift. We place the percentage of galaxies without redshift measure- 
ments at z — 00 and re-calculate Zm(mag), taking our final me- 
dian redshift estimate to be the midpoint between this upper limit 
and the measured lower limit. In the cases where the difference be- 
tween our upper and lower limit constraints are larger than the er- 
ror on the mean of the distribution a^/y/N the uncertainty on the 
median redshift is given by the upper and lower limits. For bright 
magnitudes where the redshift measurements are fairly complete 
and the number density of objects are relatively small, we place a 
statistical uncertainty on the median redshift given by / ^/N. For 
the COMBO-17 median redshift error we include the additional er- 
ror on the photometric redshift estimate where Sz/{1 + z) ~ 0.02 
for Rvega < 22 and 5z/{l + z)~ 0.05 for 22 < Rvega < 24. 

To convert the 2m (mag) from the different surveys to 
Zm(m606) we match the COMBO-17 and VVDS sources to the 
CDFS 5 epoch GOODS meos catalogue 1 Giavahsco et al. 2004). 
Note that matching the comparatively shallow redshift catalogues 
with deep 5 epoch GOODS data ensures that we are not subject to 
incompleteness in the meoe data. We then calculate the mean meoe 
in each Iab and Rvoga magnitude bin, assigning an uncertainty 



in the meoe magnitude of each bin given by aeoe/V^- Figure |8| 
shows the combined results from both surveys which are in very 
good agreement, and the best linear fit; 



-3.132 + 0.164 msoe (21.8 < meoe < 24.4). 



(16) 



To estimate the median redshift of our galaxy sample fainter than 
meoe ~ 24.4 we extrapolate the above relationship. This is justi- 
fied by the 2™, : meoe relationship measured in the Hubble Deep 
Field North (HDFN) jLanzetta et all 1 19961; iFemandez-Soto et alj 
Il999h shown (triangles) in Figure |8| where the median redshift 
error combines the photometric redshift error 52/(1 + 2) ~ 0.1 
( Fernandez- Soto et al. 1999) and the statistical uncertainty in each 
bin. Note that only the meoe > 24 points are shown for clarity, 
but there is also good agreement between COMBO-17, VVDS and 
HDFN at brighter magnitudes 

We estimate the redshift distribution for GEMS 0(2) through 



M 

= J2 N{i)n{z, meoeii))/ ^ N{i) 



(17) 



where we bin the GEMS sources into i = 1,M magnitude 
slices of mean magnitude mao(i{i), where each bin contains N{i) 
galaxies. 71(2, meoe («)) is calculated through equation J15> with 
-2o(m606) = 2m(meoe)/l-412 as estimated from equation <16> . 
The calculated 0(2) is very similar to a magnitude independent 
71(2) equation <15> . with Zm = 1-0 and as such, for simplic- 
ity when deriving the weak lensing theoretical models in the fol- 
lowing analyses, we will use a magnitude independent n{z) with 
Zm = 1.0 ± 0.1, where the error given derives from the accuracy 
of the 2m (meoe) fit, shown in Figure|8| 



7 ANALYSIS: 2-POINT STATISTICS 

In this section we use GEMS F606W data to measure the shear 
correlation function, the shear variance statistic and the shear power 
spectrum performing several diagnostic tests for systematic errors. 
We also determine an additional sampling error in order to account 
for the fact that GEMS images only one field. 



7.1 Jackknife Method 

In the analysis that follows we will o ften make use of the ja ckknife 
statistical method (see for example I Wall & Jenkin3l2003^ to ob- 
tain correlation functions with a robust estimate of the covariance 
matrix. The algorithm is quite simple, if a little time consuming. 
We are interested in the two-point correlation function which we 
first calculate from the whole survey and write as a data vector 

C — C{9i, 02, ). We then divide our sample into TV separate 

sub-regions on the sky of approximately equal area and calculate 

the correlation function Ci = Ci(6\,62, ), I ~ 1..N times 

omitting one sub-region in each calculation. Note that for a tradi- 
tional jackknife we would perform the measurement A'^ — A'^gaiaxy 
times removing a single galaxy each time, but this is computation- 
ally prohibited and provided A'^ is larger than the number of angula r 
bins, this modified jackknife method is valid iScranton et alj20o3) . 
Defining 



Cr =NC-iN- l)Ci , 



(18) 



the jackknifed estimate of the correlation function, C, is then given 
by the average C = {C*}. The jackknife estimated statistical co- 
variance of the correlation function C{9i) in angular bin i and the 



10 Heymans et al. 



correlation function C{9j) in angular bin j is given by 

{AC{e,)AC(e,)) ^ 



N{N -I) 

l=N 



(19) 



7.2 Tests for PSF contamination 

In this section we perform several diagnostic tests to determine 
whether residual PSF-related correlations remain after the PSF cor- 
rection of Section|4] Measuring the magnitude of PSF contamina- 
tion as a function of angular scale enables us to determine which 
angular scale shear correlations are free from sources of systematic 
errors and are therefore useful for cosmological parameter estima- 
tion. 



7.2.1 Star- galaxy cross correlation 

iBacon et alJi2003l) show that residual PSF-related distortions add a 
component C^^^ to the measured correlation function (7^7^) where 



c: 



(20) 



This method to test for residual PSF contamination in the data is 
similar to the cell averaged test described in Section|4]but in this 
case we look for correlations as a function of angular scale, not 
chip position, thereby revealing any PSF time variation effects that 
remain after our semi-time-dependent PSF correction has been ap- 
plied. We estimate C^J^ and associated errors using the modified 
jackknife method detailed in Section 17^1 with e* determined at 
rg = 5.9, the median galaxy size in the survey. We use A'^ = 25 
sub-regions in the jackknife estimate of M — 14 logarithmic an- 
gular bins from 6 — 0.2 arcmin to S = 25 arcmin. Figure|9|shows 
the resulting star-galaxy cross correlation functions C'l^^ and C^r^, 
compared to theoretical galaxy-galaxy shear correlation functions 
with Q,rn = 0.3, and erg = 0.75. We find that the star-galaxy 
cross correlation is consistent with zero indicating that the mea- 
surement of galaxy-galaxy shear correlations from the GEMS data 
will be free from major sources of systematics. For comparison we 
also measure the star-galaxy cross correlation when we have not 
included a correction for the distorting PSF. This reveals a corre- 
lation signal (shown dashed) that exceeds the weak cosmological 
signal that we wish to measure, stressing the importance for a good 
understanding of the PSF. 



7.2.2 E/B mode decomposition of shear correlations 

An alternative diagnostic to determine the level of systematic er- 
rors is to decompose the shear correlation function into 'E-modes' 
and 'B-modes' ( Crittenden et al. 2002). Weak gravitational lensing 
produces gradient curl-free distortions (E-mode), and contributes 
only to the curl distortions (B-mode) at small angular scales, 6 < 1 
arcmin, due to source redshift clustering ( Schneider et al. 2002). A 
significant detection of a B-type signal in weak lensing surveys is 
therefore an indication that ellipticity correlations exist from resid- 
ual system atics within the data a nd/or from intrinsic galaxy align- 
ments fsee lHevmans et all i2004l for a discussion of the latter). 

Defining the sum and difference of the tangential and rotated 
correlation functions. 
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Figure 9. Star-galaxy cross correlation functions C'^^^ (circles) and Crr^ 
(crosses) compared to theoretical galaxy-galaxy shear congelation func- 
tions (7r7r)9 (upper curve) and (7474)9 (lower curve) with = 1.0, 
f^m = 0.3, and trg = 0.7. For comparison we also measure the star- 
galaxy cross con'elation in the absence of PSF corrections, shown dashed, 
where for clarity we plot ^ (Cj^° -I- C^r'^ ) . 



aW = (7t7t)fl±(7r7r)e, (21) 

ICrittenden et alj i2002ft show that the shear correlation functions 
can be decomposed into the following E- and B-type correlators. 



where 



- 120^ 



(22) 



e-W- (23) 



Our data extends at maximum to 9 = 28 arcmin necessitating 
the use of a fiducial cosmological model to complete the integral. 
This prevents cosmological parameter estimation directly from the 
E-modes, but as the variation in the sum of the model depen- 
dent part of the integral is small ~ 10~^, this method is still a 
valid diagnostic test for residual systematics within the dat a. Note 
that the mass aperture statistic E/B decomposition ( Schnei der et alj 
ll998) does not suffer from this limitation as the integral over 
the shear correlation function spans from — > 29. This integral 
range does however introduce its own problems (see|Massev et ^ 
2004 for a discussion) and limits the analysis to small scale power 
( Van Waerbeke et all2005h . We therefore choose to use the E and 
B type correlators pur ely to test for B-type system atic errors within 
our data, although see IVan Waerbeke et alJ<2005l) for a method that 
uses the mass aperture statistic to calibrate (9) for cosmological 

parameter esti mation. 

Following iPen et alj i2002h we define the 2n element vector 
= (^+(0), ^_ (6)) which we compute from GEMS, binning the 
data finely into n — 2000 intervals of 0.9 arcsec (equivalent to a 
separation of ~ 30 ACS pixels). The E/B correlators are then given 
by a 2?! element vector 



(24) 



where T is a 2n^ transformation matrix defined by equation <23> . 
To reduce noise we re-bin ^f*^ into M — 14 logarithmic intervals 
from 9 = 0.2 arcmin to 6* = 25.0, represented by an Al x 2n 
projection operator ^J?'" = Pki^f^. 

To calculate errors on we first calculate the covariance 
matrix of the raw correlators Lij = (A^iA^j) using the jackknife 
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Figure 10. E/B mode decomposition of the stiear correlation function mea- 
sured from tlie GEMS observations. Tlie upper panel shows g^'" and the 
fiducial ACDM theoretical £,^(9) model where erg = 0.7 and the me- 
dian source redshift Zm = 1.0. The lower panel shows 5^'" that is consis- 
tent with zero on all angular scales and can be compared to the theoretical 
model ^^{0) shown dotted. Using a different fiducial ACDM cosmologi- 
cal model to calculate has a small effect (at the level of 10~^) which 
can be seen from the dotted curves where we have assumed erg = 0.6 
(lower curve E-mode, upper curve B-mode) and eg = 1.0 (upper curve 
E-mode, lower curve B-mode) 



metho d detailed in Section ITTI We find, in contrast to |Pen et'all 
i2002h . that our slightly wider correlation bins (A9 = 0.9 arcsec 
compared to A9 = 0.6 arcsec) are correlated. Note that the modi- 
fied jackknife method in principle requires N > n sub-regions to 
calculate Lij, but as we re-bin Lij into M = 14 broader angu- 
lar scales to estimate the final errors on our jackknife method, 
which is computationally time limited to A'^ = 100, is still valid. 
The binned cov ariance matrix of the E/B correlators is given by 
jPen et all20ol . 



(25) 



/ r bin 
^kk ■ 



where the errors on are given by 

We investigate E- and B-mode correlations in the GEMS and 
GOODS data separately with the measurement from the GEMS 
data shown in Figure^] For this E/B mode decomposition we have 
used a fiducial ACDM cosmological model to complete the inte- 
gral of equation <23> with as — 0.7 and with a median source 



redshift for our galaxies of Zm = 1.0. Using a different fidu- 
cial ACDM cosmological model has a small effect (at the level 
of 10^^) which can be seen from the dotted curves in Figure [Tol 
where we have assumed ag = 0.6 and ag — 1.0. We find that 
the E-modes are in good agreement with the fiducial cosmological 
model and that the B -modes for GEMS are consistent with zero on 
all scales 6 > 0.2 arcmin. Whilst finding this result very encour- 
aging we note that our B-modes measured at small angular scales 
are very strongly correlated and noisy such that in the worse-case 
scenario, represented by the upper end of the error bars shown, the 
B-modes exceed the signal that we wish to measure. This motivates 
our desicion to conservatively limit our shear correlation function 
analysis to angular scales 6 > 0.65 arcmins within the GEMS data 
where we can be confident that the signal we measure is cosmolog- 
ical and not systematic. 

The PSF model for the GOODS data is determined from ap- 
proximately half the number of stars which were used to deter- 
mine the two semi-time-dependent GEMS PSF models, as GOODS 
spans approximately one quarter of the area of the GEMS observa- 
tions. We would therefore expect to find a poorer PSF correction 
with the GOODS data which is seen with the presence of non-zero 
B-modes at angular scales 6^ < 1 arcmin. We therefore only include 
GOODS data in our shear correlation function analysis for angular 
scales 9 > 1 arcmin. 



7.3 The shear correlation function 

Having shown in the previous section that we are not contaminated 
by significant non-lensing correlations, we can now measure the 
weak lensing shear correlation function {"/{0)-y{9 + A6)). To ob- 
tain results that are independent of the initial frame of reference 
we measure the tangential and rotated shear correlation functions, 
equation ^3} and equation ^4) respectively. These can be estimated 
from the data by 



Npairs ^-^ 

pairs 



7'(x) 7t(x + 6»), 



(26) 



where the tangential shear 74 and rotated shear jr are given by 
equation 74 — 7™* and 7r = 72°* where the rotation angle (f) 
is now defined to be the angle between the x axis and the line join- 
ing each galaxy pair. Note this rotation follows the initial rotation 
that sets the full GEMS shear catalogue into the same reference 
frame. We also calculate the cross-correlation function E['yt'yr]e, 
which the parity invariance of weak lensing predicts to be zero. 

We calculate the shear correlations using the modified jack- 
knife method detailed in Section lTTI We use TV = 25 sub-regions 
in the jackknife estimate of 71/ = 11 logarithmic angular bins from 
6 = 0.65 arcmin to 61 = 25 arcmin and we include data from 
the GOODS area only for angular scales 9 > 1 arcmin. This en- 
sures that the shear correlation measurement at small angular scales 
^ < 1 arcmin is not contaminated by the small-scale non-lensing 
distortions found in the GOODS data. We show the resulting jack- 
knife estimates for the tangential and rotated shear correlation func- 
tions in Figure^3 The theoretical AC DAI models over-plotted are 
calculated from equation ^3} and equation J4j with a median galaxy 
redshift 2™ = 1.0 and ag = (0.6, 0.7, 1.0). Note that we find the 
cross-correlation i5[7t7r]e to be consistent with zero on all scales, 
as expected, supporting our findings that we are not contaminated 
by significant non-lensing distortions. 
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Figure 11. Shear correlation functions -E[7t7t]e (upper panel) and 
E[^r'yr\g (lower panel) estimated from the GEMS data using a modified 
jackJoiife technique. The GOODS section of the GEMS survey is only in- 
cluded in measurements of i?[77]a for > 1.0 arcmin so as not to include 
the non-lensing B-mode systematics seen in the GOODS data at smaller an- 
gular scales. Over-plotted are theoretical ACDM models {-ft'yt)e (upper 
panel) and (7r7r)e (lower panel) with erg = 0.6 (dashed lower), = 0.7 
(solid) and ag = 1.0 (dashed upper). 



7.3.1 Sampling variance 

The GEMS mosaic samples only one area of sky and as such our re- 
sults are subject to additional sampling variance errors. To address 
this issue we have created 100, 28 x 28 arcmin, Gaussian realisa- 
tions of a shear power spectrum calculated for aag = 0.75, ACDM 
cosmology for sources with Zm = i. We measure the shear correla- 
tion function from each realisation, populating the shear field with 
60 circular galaxies per square arcmin. The variance we measure 
between the results from each realisation then provides us with an 
estimate of the sampling variance. We calculate the sampling co- 
variance matrix from the 100 independent realisations and add this 
to the jackknifed covariance matrix measured from the GEMS data, 
slightly overestimating the sampling variance error on small scales. 
The reader should note that this method is somewhat cosmology 
dependent but it is sufficient to estimate the amplitude of sampling 
variance for our cosmological parameter constraints in section [8| 
Future work will investigate the impact of non-Gaussianity on our 
sampling variance est imation using the cosmological ray-tracing 
N-body simulations of IWhitd i200g 



7.4 The shear variance statistic 

Space-based lensing surveys to date and early ground-based lens- 
ing results focused on the shear variance statistic, equation jsj, to 
analyse the data. This statistic produces the highest signal-to-noise 
measurement of weak lensing shear and can be estimated from the 
data by splitting the sample into N circular cells of radius 9 and cal- 



culating the shear variance in excess of noise (see lBrown et alj2003l 
for a minimum variance estimator). As discussed in Section l7.2.2l 
measuring the B-mode of the shear correlation function allows one 
to select angular scales above which one can be confident that the 
presence of non-lensing distortions are insignificant. Very small 
scale systematic distortions are poorly understood and are sucess- 
fuUy ignored by the shear correlation statistic. For the shear vari- 
ance statistic however, small scale non-lensing distortions are in- 
cluded in the measurement at all angular scales biasing the shear 
variance statistic. 

To assess the impact of our residual small scale non-lensing 
distortions on the shear variance statistic we can, in a similar fash- 
ion to Section l7.2.2l decompose the signal into its E-mode and B- 
mode components. Sc hneider et al. t,200Z) show that the shear vari- 
ance of the E- and B-mode can be obtained in terms of the shear 
correlation functions (^+, equation <2U . through 



/ 2\ E 

(7 )e 



(27) 



/ 2\ B 

(7 )e = 



26(2 



(28) 



where S+ and S- are given in equation (39) and equation (42) of 
ISchneider eTal] j2002l) . Note that 5*- does not cut off at finite sep- 
aration and as such one needs to include a fiducial cosmological 
model to complete the integral, as in Section ^L22i We calculate 
{7^)f and {'y^)f following the method of Pen et al. (2002), de- 
tailed in Section l7.2.2l where the transformation matrix T of equa- 
tion <24> is now defined by equation <27t and equation <28t . 

Fi sure [T2I shows the result of our E/B mode decomposition of 
the shear variance statistic revealing significant B-modes on scales 
6 < 1.5 arcmin. Replacing the small scale 5+ and ^_ measure- 
ments (0 < 0.2 arcmin) with a theoretical model, we find B-modes 
consistent with zero on all scales. The B-modes shown in Fieure fT^ 
therefore result from very small scale strong non-lensing distor- 
tions that bias the shear variance statistic even at larger angular 
scales. Non-lensing distortions are likely to contribute equally to 
the measured E-mode making the signal appear to favour a higher 
value for as, when compared to the large scale B-mode free shear 
variance measurements. These large angular scales have previously 
been unmeasureable from space-based data. This demonstrates that 
in order to obtain reliable cosmological parameter constraints it is 
vital to perform an E/B mode decomposition of the statistic that 
will be used to constrain cosmology. We will not use the shear 
variance statistic to constrain cosmological parameters, even on B- 
mode free scales, preferring the shear correlation function statistic 
in Section rOl and the shear power spectrum that we determine in 
the following section. 



7.5 Shear power spectrum 

In addition to the shear correlation function and shear variance 
statistic of the previous sections, we also quantify the two-point 
statistics of the shear field by directly measuring its power spec- 
trum, Pk, equation 0. Power spectrum estimation from cosmo- 
logical data sets is a well-studied problem in the context of mea- 
suring the statistical properties of the CMB (see lEfstathioijE004l 
for an overview) and the methods developed in this field are com- 
pletely applicable to measuring power spectra from weak lensing 
data sets . Here, we use a maximum likelihood estimator (see for 
example jBond et ai]|l998h to reconstruct the power spectrum of 
the shear field observed in the GEMS data. Our approach is based 
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Figure 12. E/B mode decomposition of tlie shear vaiiance. Tlie upper panel 
shows (7^)!' and the tiducial ACDM theoretical {■^'^) model (solid) where 
CTg = 0.7 and the median source redshift = 1-0. The lower panel 
shows {'y^)f that is consistent with zero on angular scales 9 > 1.5 arcmin. 
Using a different fiducial ACDM cosmological model to calculate (7^)|'^ 
has a small effect (at the level of 3 X 10^^) which can be seen from the 
dotted curves where we have assumed ag = 0.6 (lower curve E-mode, 
upper curve B-mode) and ag = 1.0 (upper curve E-mode, lower curve B- 
mode). The significant non-lensing B-modes al 9 < 1.5 ai'cmin indicate 
residual small scale systematics that, with the shear variance statistic, ai'e 
included at all angular scales. These non-lensing distortions also contribute 
to the measured E-mode making the signal appear to favour a high value 
for erg (upper panel: ACDM erg = 1.0 theoretical model over-plotted dot- 
dashed). 



tion l7.2.2l and the E/B shear variance measurements of Section l7!4l 
where a fiducial cosmological model is needed to complete inte- 
grals over the infinite correlation function thereby invalidating their 
sole use for cosmological parameter estimation. 

To apply the estimator to the GEMS data, we bin the galaxy 
shear distribution into 30 x 30 equal-size pixels of ~ 1 square 
arcmin. Writing this pixelised shear distribution as a vector, d, we 
then maximise the likelihood function, 



21nL(C|d) =d*C"M + Tr[InC] 



(29) 



using a Newton-Raphson scheme, as a function of the band-powers 
of the three power spectra, P"", P'^'^ and P"'^ . Here, C is the 
data covariance matrix which is the s um of the cosmological signal 
(equation (21) of iBrown et alJ ^2bo3^ ^ and a noise term. 



N : 



iVoix 



I, 



(30) 



where "/nus and Npi^ are the root mean square shear and occu- 
pation number within each pixel respectively. The errors and co- 
variance of our final band-powers are approximated as the inverse 
Fisher matrix, which is an excellent approximation provided that 
the likelihood function is sufficiently Gaussian in the band-powers. 

Figure fT3l shows the results of applying the maximum like- 
lihood estimator to the GEMS data along with a theoretical shear 
power spectrum for a ACDM model with Qrn = 0.3 and as = 0.8 
with which we find reasonable agreement. The measurements of 
the B-mode spectrum are mostly consistent with zero although 
there is a significant detection of E-B cross-correlations on medium 
scales. We suspect that these come from the time-variation of the 
PSF that we have only partially accounted for with our semi-time 
dependent PSF modelling. Our PSF models are designed to de- 
crease the average stellar ellipticity to zero and therefore when 
averaging over the whole survey, as in the measurement of the 
shear correlation function, the residual PSF contamination is zero, 
as can be seen in Figure |9| For the shear power spectrum mea- 
surement however, where the field is decomposed into its Fourier 
components, the time-dependent PSF contamination can be identi- 
fied. Note that from the covariance of the P"" measurements we 
find that our band-power measurements are almost independent of 
one another, apart from the slight anti-correlation of neighbouring 
bands which is a natural consequence of the maximum likelihood 
estimator. 



on the prescription of IHu & White! i200ll) who proposed recon- 
structing the three power spectra, P"^"^, pl^^ and P"**^ as a series 
of step-wise 'band-powers' where the quantity £{£ + l)P'-'/27r is 
approximated as a constant within each band. P*^*^ is the power 
spectrum of the B-modes while P'''' is the cross power-spectrum 
between the E- and B-modes. The maximum likelihood method 
automatically accounts for irregular survey geometries, pixeliza- 
tion effects and produces error estimates, via the Fisher Informa- 
tion matrix (see for example Teamarket al. 1997), which include 
sampling variance and shot noise. .Hu & White (.2001.) have tested 
the maximum likelihoo d estimator on both Gaussian and N-body 
simulations, while Brown et alJ l200'3 ) have tested it on Gaussian 
simulations on scales similar to the GEMS data and have applied 
the estimator to the COMBO- 17 weak lensing data set. 

The maximum-likelihood decomposition only necessitates the 
use of a fiducial cosmological model to estimate the significance 
of the result. This is in comparison to the E/B correlators of Sec- 



8 COSMOLOGICAL PARAMETER ESTIMATION 

Having measured the 2-point statistics of the shear field within 
GEMS, we can now compare these measurements with theoret- 
ical predictions in order to place joint constraints on the matter 
density of the Universe Q,rn and the normalisation of the matter 
power spectrum ag. We do this using both our correlation func- 
tion measurements from Section FOl and the power spectrum esti- 
mates from Section 1731 We use equation Q to calculate our the- 
oretical shear power spectra and equation ^3} and equation ^4} to 
calculate our theoretical correlation functions for a variety of cos- 
mological mod els. For these calculatio ns we have used the trans- 
fer function of 'Eisens tein & Hul ^1999") for the dark matter power 
spectrum with an initial power spectruin slope of n = 1. To pro- 
duce the no n-linear power spec trum from this, we use the fitting 
formulae of ISmith et alj i2003h and we fix Qm + = 1. We 
also use the form of equation <15> for the input redshift distribu- 
tion of source galaxies. We consider models in the following re- 
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Figure 13. The cosmic shear power spectra from GEMS. Plotted on a linear- 
log scale are P'^"' (circles), P'''' (crosses) and P'^^ (triangles) in five band- 
averaged band-powers as a function of multipole, I. The errors bars are es- 
timated from the Fisher matrix and the Pl^^^ and P"'' have been slightly 
horizontally displaced for clarity. The solid curve is the shear power spec- 
trum estimated for a erg = 0.8 normalised ACDM model. 



where x = x(ct8, ^m,Ho, Zm) are the theoretical correlation func- 
tions ordered in a similar manner to our data vector. V = (dd^) is 
the sum of the covariance matrix of our correlation function mea- 
surements as estimated from the data using equation <19t and a 
sampling covariance matrix as detailed in section rf.B.ll The fitting 
of the power spectrum measurements is done in a similar fashion 
where V, the covariance matrix of the band-power measurements, 
comes from a Fisher error analysis that automatically includes sam- 
pling variance. After calculating values for each of our theoret- 
ical models, we marginalise over the Hubble constant. Ho with a 
prior set by the I'" year WMAP results (Ho = 72 ± 5 km s~^ 
Mpc"^ Soeraeletal. 1 2003)). We also marginalise over the me- 
dian redshift of the source galaxies, Zm using Zm ~ 1.0 ± 0.1 as 
estimated in Section|5| 

The resulting constraints in the — Qm plane for both the 
shear correlation function and shear power spectrum estimators are 
shown in Figure 1141 We find good agreement between the con- 
straints obtained using the two different measures: for the corre- 
lation function measurements, we find 



as{nm/0-3)°-''^ = 0.68 ± 0.13 



(33) 



while using the power spectrum analysis, we find a slightly higher 
value of 



asinm/O.ay-^^ = 0.72 ± 0.11. 



(34) 
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Figure 14. The likelihood surface of ag and Qm from GEMS as calcu- 
lated using the shear correlation function (dark contour) and using the shear 
power spectrum (light contours). We plot the Icr and 2cr confidence regions 
for the shear power spectrum, and, for clarity, the Icr confidence regions for 
the shear correlation function 



9 CONCLUSION 

In this paper we have presented the detection of weak gravitational 
lensing by large-scale structure in the GEMS survey, demonstrating 
that our shear correlation signal is uncontaminated by significant 
non-lensing shear distortions. GEMS, imaged by the ACS on HST, 
spanning 795 square arcmin, is the largest contiguous space-based 
mosaic that has undergone a cosmic shear analysis to date. This has 
enabled us to measure cosmic shear over a large dynamic range of 
angular scales; from the small scales (6 = 0.65 arcmin) that are 
difficult to probe with ground-based surveys, up to the large scales 
(9 = 21.0 arcmin) that were previously inaccessible to space-based 
surveys. Our careful analysis, where we have considered forms of 
selection bias, centroid bias and calibration bias, geometric shear 
distortions and PSF contamination, has yielded an unbiased^ mea- 
surement of the shear correlation function uncontaminated by non- 
lensing 'B-mode' distortions. This has allowed us to set joint con- 
straints uncontaminated by major sources of systematic errors on 
the matter density of the Universe flm and the normalisation of the 
matter power spectrum ag finding ag = 0.73 ± 0.13 for WMAP 
constrained D,m = 0.27 fSDergel et al]|2003). It is interesting to 
note that the GEMS cosmological parameter c onstraints are very 
similar to those froin the COMBO- 17 survey (.Brown et alJl2003l: 
iHevmans et all2004l) . a deep multi-colour survey which spans ~ 4 



gions of parameter space: 0.3 ^ erg ^ 1.5 ; 0.1 ^ Qm ^ 1-0 ; 
64 s: Ho s; SOkms^^Mpc"^ and 0.9 < 2„ 1.1. 

Writing our correlation function measurements as a data vec- 
tor, 

{Ci{9i),...,Ci{9„),C2{9i),...,C2{9n)}, (31) 
for each theoretical model, we calculate 



Our measurement is unbiased if we assume that the KSB-l- method ap- 
plied provides us with an unbiased estimate of galaxy shear 7, which has 
been shown to be true with ground-based data jErben et al. 20013. The equal 
galaxy shear (on average) measured in our F606W and F850LP data sug- 
gests that the impact of strongly non-Gaussian space-based PSFs on the 
KSB+ method is small, supporting its use as an unbiased shear estimator 
in this paper. This will be investigated further with sheared spaced-based 
image simulations in a forthcoming paper. 



Weak lensing with GEMS 15 



times the area of GEMS. This results from the higher number den- 
sity of resolved galaxies in space-based data and the higher signal- 
to-noise measurements of galaxy shear which are achievable with 
higher resolution data (Brown et al. in prep). 

We have presented a thorough discussion on the anisotropic 
ACS PSF that, for the first time with a space-based weak lensing 
analysis, we have been able to characterise directly from our data 
without having to assume long-term PSF stability. Long-term PSF 
stability has been tested with the GEMS data and shown to be true 
for the ACS only above the ~ 5% level. With previous space-based 
analyses it has often been assumed that temporal changes in the 
mean PSF pattern can be accommodated by a low order correc- 
tion. Percent level accuracy in the ACS PSF correction would how- 
ever be difficult to achieve with this assumption as this work has 
shown that the both the PSF strength and pattern exhibits temporal 
variation. We have identified PSF temporal variation on the level 
of a few percent finding consistent behaviour between the F850LP 
imaging and F606W imaging, even though the F850LP PSF is quite 
different from the F606W PSF. We have tested the success of our 
PSF correction by measuring the star-galaxy cross-correlation and 
the B-type shear correlator which were both found to be consis- 
tent with zero on angular scales 9 > 0.2 arcmin. Our semi-time- 
dependent method for PSF modelling therefore adequately corrects 
for the varying PSF distortion when we consider weak lensing shear 
correlations as a function of relative galaxy position averaged over 
the whole GEMS mosaic. When we measure the shear power spec- 
trum however, a statistic which is dependent on galaxy shear as a 
function of absolute galaxy position, we find a significant detec- 
tion of E-B cross-correlations most likely revealing the impact of 
not producing a fully time-dependent PSF correction model. It is 
currently unclear where the variation in the GEMS ACS PSF arises 
but its presence, also seen bv Jee et al.. (.2005 ) and Rhodes et. al. 
(in prep), suggests that future HST cosmic shear surveys should be 
preferentially observed in sequence to minimise the impact of PSF 
instabilities. 

We have measured the commonly used top-hat shear variance 
statistic, performing an E/B mode decomposition. We find signif- 
icant non-lensing B-mode distortions at angular scales 6 < 1.5 
arcmin in contrast to the E/B decomposition of the shear corre- 
lation function where the B-modes were found to be consistent 
with zero at angular scales 6 > 0.2 arcmin. This shows that 
the top-hat shear variance statistic becomes contaminated by very 
small scale non-lensing distortions out to fairly high angular scales, 
strongly biasing the final result. Note that this effect is also seen in 
[Van Waerbeke et al. ( 2005). The shear correlation function does not 
suffer from this contamination as the very small scale non-lensing 
correlations are removed from the analysis and it is therefore this 
statistic along with the shear power spectrum that we favour for 
cosmological parameter estimation. We urge future cosmic shear 
studies to perform E/B mode decompositions to test for non-lensing 
distortions and employ statistical analyses other than the easily bi- 
ased top-hat shear variance statistic. 
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Figure 15. Comparison of the total shear correlation function E[yy]g as 
measured from GEMS along with the most up-to-date shear correlation 
measurements from the other groups indicated. Over-plotted ai'e theoreti- 
cal AC DM models for a Zm = 1 survey with erg = 0.7 (lower) and 
CTg = 1.0 (upper). Note all data points and errors have been scaled to a 
Zm = 1 survey using a redshift scaling. 
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9.1 Comparison with other cosmic shear surveys 

Figures fTsl and fT6l compare the GEMS results with the most up- 
to-date results from other cosmic shear surveys that have placed 
constraints on as and Sl^ . The cosmic shear signal scales with the 
depth of the survey and so we have introduced a median redshift 



Figure 16. Comparison of the top-hat shear variance E['y'^]g as measured 
from GEMS along with the most up-to-date top-hat shear variance mea- 
surements from the other groups indicated. We show only the B-mode free 
GEMS top-hat shear variance results with 9 ^ 1.5 arcmin. Over-plotted ai'e 
theoretical AC DM models for a Zm = 1 survey with erg = 0.7 (lower) 
and erg = 1.0 (upper). Note all data points and errors have been scaled to a 
Zm, = 1 survey using a z^*^ redshift scaling. 
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scaling* of the data points and errors to bring the different results in 
line with a survey of median redshift Zm ~ 1-0. For the shear cor- 
relation function (Figure [TsJ we scale by z^, as suggested by the 
numerical simulations of Barber (2002) and fo r the top-hat shear 
variance (Figure fT6l we scale by z;^^-' (Rho des et alj|2004 . We 
preferentially show the shear correlation function for surveys that 
have measured the top-hat shear variance in addition, due to our 
concern with the use of the top-hat shear variance statistic. We only 
show the B-mode free GEMS top-hat shear variance results with 
6 ^ 1.5 arcmin. These comparisons show broad agreement be- 
tween the shear correlation measurements and a poorer agreement 
between the top-hat shear variance measurements. As discussed in 
Section IV^l the top-hat shear variance appears to be easily con- 
taminated at large scales by small scale systematic errors and we 
propose that this contamination, not always quantified, is at least 
a partial cause of the differing results. Other possibilities are po- 
tential calibration biases arising from differences betw een the vari- 
ous sh ear measureme nt methods (compare for example lErben et all 
I2OOII and lBacon etaL200 1). differences in the median redshift de- 
termination and sampling variance. 

The results shown in Figure fTsl and Figure fT6l vield measure- 
ments of (78 ranging from as — 0.7 to erg — 1.1 for a value of 
fim = 0.3. This can be compared w ith results from the WMAP 
CMB experiment ( Soergel et ai]l2003lh that finds erg = 0.9 ± 0.1 
from the WMAP data alone and as = 0.84 ± 0.04 when the 
WMAP data is combined with other data sets. Results from clus- 
ter abundance measurements range from erg = 0.7 to as = 1.0 
(see[PierDaoli et al. 2003, who find as = 0.77 ± 0.05, and refer- 
ences therein). Our measurement is at the lower end of all these 
results which we may expect in light of the fac t that the CDFS i s 
a factor of two under-dense in massive galaxies ^Wolfetal."2003|). 
If we assume that massive galaxies trace the underlying dark mat- 
ter distribution, then we would expect a low measurement of as 
from this field when compared to the global erg value. Combining 
GEMS data with other wide-field space-based mosaics, such as the 
COSMOS survey^ and the ACS pure parallel survey, will reduce 
the effects of sampling variance in order to obtain a good estimate 
of the Universal value of as from HST. 



10 ACKNOWLEDGEMENTS 

This work is based on observations taken with the NASA/ESA 
Hubble Space Telescope, which is operated by the Association of 
Universities for Research in Astronomy, Inc. (AURA) under NASA 
contract NAS5-26555. Support for the GEMS project was provided 
by NASA through grant number GO-9500 from the Space Tele- 
scope Science Institute, which is operated by AURA for NASA, 
under contract NAS5-26555. CH and HWR acknowledge financial 
support from GIF. MLB and CW were supported by PPARC fel- 
lowships. EFB and SFS acknowledge financial support provided 



* For the results from lChang et alj |2004) we scale assuming the me dian 
radio s ource redshift to be Zm = 2.0. For the results from Casertano et alj 
l2003ft we convert the measurement from the top-hat variance in square 
cells to the top-h at variance in circular cells using a 1/\/tt scaling 
iBaco n et al.'200(j)^ and then scale the results using the median redshift de- 
rived by Refregi er et alj l2002h for the same data set. For the ljarvis et alj 
12003) results we scale by Zm = 0.6 (private communication) and for all 
other results, we use the quoted median redshift. 
5 www. astro. caltech.edu/^ cosmos 



through the European Community's Human Potential Program un- 
der contract HPRN-CT-2002-00316, SISCO (EFB) and HPRN-CT- 
2002-00305, Euro3D RTN (SFS). SJ and DHM acknowledge sup- 
port from NASA under LTSA Grant NAG5-13063(SJ) and under 
LTSA Grant NAG5-13102(DHM), issued through the Office of 
Space Science. We thank Justin Albert, Colin Cox, Richard Ellis, 
Yannick Mellier, John Peacock, Jason Rhodes, Tim Schrabback, 
Ludovic Van Waerbeke and Martin White for many useful discus- 
sions, and also the referee for his/her helpful comments. 



REFERENCES 

Bacon D., Massey R., Refregier A., EUis R., 2003, MNRAS, 344, 
673 

Bacon D., Refregier A., Elhs R., 2000, MNRAS, 318, 625 
Bacon D. J., Refregier A., Clowe D., Ellis R. S., 2001, MNRAS, 

325, 1065 
Barber A. J., 2002, MNRAS, 335, 909 

Barden M., Rix H.-W., Somerville R. S., Bell E. F, HauBler B., 
Beckwith S. V. W., Borch A., Caldwell J. A. R., Heymans C, 
Jahrtke K., Jogee S., Mcintosh D. H., Meisenheimer K., Peng 
C. Y, Sanchez S. F, Wisotzki L., Wolf C, 2004, ApJ submitted, 
astroph/0502416 

Bartelmann M., Schneider P., 2001, Physics Reports, 340, 291 

Baugh C. M., Efstathiou G., 1994, MNRAS, 267, 323 

Bernstein G. M., Jarvis M., 2002, AJ, 123, 583 

Bertin E., Arnouts S., 1996, A&AS, 117, 393 

Bond J. R., Jaffe A. H., Knox L., 1998, Phys. Rev. D, 57, 2117 

Brown M., Taylor A., Bacon D., Gray M., Dye S., Meisenheimer 
K., Wolf C, 2003, MNRAS, 341, 100 

Casertano S., Ratnatunga K. U., Griffiths R. E., 2003, ApJL, 598, 
L71 

Chang T, Refregier A., Helfand D. J., 2004, ApJ, 617, 794 
Contaldi C, Hoekstra H., Lewis A., 2003, Phys. Rev. Lett., 90, 
303 

Cox C, Gilliland R. L., 2002, in The 2002 HST Calibration Work- 
shop The Effect of Velocity Aberration on ACS Image Process- 
ing, pp 58— F 

Crittenden R., Natarajan R., Pen U., Theuns T, 2002, ApJ, 568, 
20 

Efstathiou G., 2004, MNRAS, 349, 603 
Eisenstein D. J., Hu W., 1999, ApJ, 511, 5 

Erben T, Van Waerbeke L., Bertin E., Mellier Y, Schneider P., 
2001, A&A, 366,717 

Fernandez-Soto A., Lanzetta K. M., Yahil A., 1999, ApJ, 513, 34 

Ford et al. H. C, 2003, in Proceedings of the SPIE, Volume 4854, 
pp. 81-94 (2003). Overview of the Advanced Camera for Sur- 
veys on-orbit performance, pp 81-94 

Giavalisco et al. M., 2004, ApJL, 600, L93 

Hammerle H., Miralles J.-M., Schneider P., Erben T, Fosbury 
R. A. E., Freudling W., Pirzkal N., Jain B., White S. D. M., 2002, 
A&A, 385, 743 

Hamana T, Miyazaki S., Shimasaku K., Furusawa H., Doi M., 
Hamabe M., Imi K., Kimura M., Komiyama Y, Nakata F, Okada 
N., Okamura S., Ouchi M., Sekiguchi M., Yagi M., Yasuda N., 
2003, ApJ, 597, 98 

Heymans C, Brown M., Heavens A., Meisenheimer K., Taylor 
A., Wolf C, 2004, MNRAS, 347, 895 

Heymans C, Heavens A., 2003, MNRAS, 339, 711 

Hirata C, Seljak U., 2003, MNRAS, 343, 459 



Weak lensing with GEMS 17 



Hoekstra H., Franx M., Kuijken K., Squires G., 1998, ApJ, 504, 
636 

Hoekstra H., Yee H., Gladders M., 2002, ApJ, 577, 595 
Hu W., White M., 2001, ApJ, 554, 67 

Jain B., Taylor A., 2003, Physical Review Letters, 91, 141302 
Jarvis M., Bernstein B. J. G., Dolney D., 2005, ApJ in press, as- 
troph/0502243 

Jarvis M., Bernstein G., Jain B., Fischer P., Smith D., Tyson J., 

Wittman D., 2003, ApJ, 125, 1014 
Jee M. J., White R. L., Bemtez N., Ford H. C., Blakeslee J. P., 

Rosati R, Demarco R., lUingworth G. D., 2005, ApJ, 618, 46 
Kaiser N., 2000, ApJ, 537, 555 

Kaiser N., Squires G., Broadhurst T, 1995, ApJ, 449, 460 

King L., Schneider R, 2002, A&A, 396, 411 

Koekemoer A. M., Fruchter A. S., Hack W. J., Hook R. N., 2003, 

in HST Calibration Workshop (STScI: Baltimore) Multi-drizzle. 

p. 337 

Krist J., 2000, TinyTim User's Guide v4. (STScI) 
Lanzetta K. M., Yahil A., Fernandez-Soto A., 1996, Nature, 381, 
759 

Le Fevre et al. O., 2004, A&A, 428, 1043 

Luppino G. A., Kaiser N., 1997, ApJ, 475, 20 

MacDonald E. C., Allen P., Dalton G., Moustakas L. A., Heymans 
C, Edmondson E., Blake C., Clewley L., Hammell M. C, Old- 
ing E., Miller L., Rawlings S., Wall J., Wegner G., Wolf C., 2004, 
MNRAS, 352, 1255 

Maoli R., Van Waerbeke L., Mellier Y., Schneider P., Jain B., 
Bernardeau F, Erben T, 2001, A&A, 368, 766 

Massey R., Bacon D., Refregier A., Ellis R., 2004, MNRAS sub- 
mitted, astroph/0404195 

Massey R., Refregier A., 2004, MNRAS submitted, as- 
troph/0408445 

Massey R., Rhodes J., Refregier A., Albert J., Bacon D., Bernstein 
G., EUis R., Jain B., McKay T., Perlmutter S., Taylor A., 2004, 
AJ, 127, 3089 

Meurer G. R., Lindler D. J., Blakeslee J., Cox C. R., Martel A., 
Tran H. D., Bouwens R., Ford H. C, Clampin M., Hartig G. F, 
Sirianni M., De Marchi G., 2003, in Proceedings of the SPIE, 
Volume 4854, pp. 507-514 (2003). CaUbration of geometric dis- 
tortion in the ACS detectors, pp 507-514 
Miralles J.-M., Erben T, Hammerle H., Schneider P., Freudling 

W, Pirzkal N., Fosbury R. A. E., 2005, A&A, 432, 797 
Park Y, Casertano S., Ferguson H. C, 2004, ApJL, 600, L159 
Pen U.-L., Van Waerbeke L., Mellier Y, 2002, ApJ, 567, 31 
Peng C. Y, Ho L. C, Impey C. D., Rix H.-W., 2002, AJ, 124, 266 
PierpaoU E., Borgani S., Scott D., White M., 2003, MNRAS, 342, 
163 

Refregier A., Bacon D., 2003, MNRAS, 338, 48 

Refregier A., Massey R., Rhodes J., Ellis R., Albert J., Bacon D., 

Bernstein G., McKay T, Perlmutter S., 2004, AJ, 127, 3102 
Refregier A., Rhodes J., Groth E. J., 2002, ApJL, 572, L131 
Rhodes J., Refregier A., CoUins N. R., Gardner J. P., Groth E. J., 

HiUR.S., 2004, ApJ, 605,29 
Rhodes J., Refregier A., Groth E. J., 2000, ApJ, 536, 79 
Rhodes J., Refregier A., Groth E. J., 2001, ApJL, 552, L85 
Rhodes J., Refregier A., Massey R., Albert J., Bacon D., Bernstein 

G., EUis R., Jain B., Kim A., Lampton M., McKay T., SNAP 

Collaboration 2004, Astroparticle Physics, 20, 377 
Riess A., 2002, in The 2002 HST Calibration Workshop Growth 

of Hot Pixels and Degradation of CTE for ACS. pp 47-+ 
Riess et al. A. G., 2004, ApJL, 600, L163 
Rix H.-W., Barden M., Beckwith S. V. W., Bell E. F., Borch A., 



Caldwell J. A. R., HauBler B., Jahnke K., Jogee S., Mcintosh 
D. H., Meisenheimer K., Peng C. Y, Sanchez S. R, Somerville 
R. S., Wisotzki L., Wolf C, 2004, ApJS, 152, 163 
Schneider R, Van Waerbeke L., Jain B., Kruse G., 1998, MNRAS, 

296, 873 

Schneider P, Van Waerbeke L., MeUier Y, 2002, A&A, 389, 741 
Scranton et al. R., 2002, ApJ, 579, 48 

Smith R., Peacock J., Jenkins A., White S., Frenk C, Pearce R, 
Thomas P., Efstathiou G., Couchman. H., 2003, MNRAS, 341, 
1311 

Spergel D., Verde L., Peiris H., WMAP Collaboration 2003, 
ApJS, 148, 175 

Tegmark M., Taylor A. N., Heavens A. R, 1997, ApJ, 480, 22 
Tereno I., Dore O., van Waerbeke L., Mellier Y, 2005, A&A, 429, 
383 

Van Waerbeke L., Mellier Y, Hoekstra H., 2005, A&A, 429, 75 
Van Waerbeke L., Mellier Y, Radovich M., Bertin E., Dantel-Port 

M., McCracken H., Fevre O. L., Foucaud S., Cuillandre J., Erben 
T., Jain B., Schneider R, Bernardeau F, Fort B., 2001, A&A, 
374, 757 

Wall J. v., Jenkins C. R., 2003, Practical statistics for astronomers. 
Princeton Series in Astrophysics 

White M., 2005, Aph in press astroph/0502003 

Wolf C, Meisenheimer K., Kleinheinrich M., Borch A., Dye S., 
Gray M., Wisotzki L., Bell E., Rix H.-W, Hasinger A. C. G., 
Szokoly G., 2004, A&A, 421, 913 

Wolf C, Meisenheimer K., Rix H.-W, Borch A., Dye S., Klein- 
heinrich M., 2003, A&A, 401, 73 



